Array-representation integration factor method for high-dimensional systems
نویسندگان
چکیده
High order spatial derivatives and stiff reactions often introduce severe temporal stability constraints on the time step in numerical methods. Implicit integration method (IIF) method, which treats diffusion exactly and reaction implicitly, provides excellent stability properties with good efficiency by decoupling the treatment of reactions and diffusions. One major challenge for IIF is storage and calculation of the potential dense exponential matrices of the sparse discretization matrices resulted from the linear differential operators. Motivated by a compact representation for IIF (cIIF) for Laplacian operators in two and three dimensions, we introduce an array-representation technique for efficient handling of exponential matrices from a general linear differential operator that may include cross-derivatives and non-constant diffusion coefficients. In this approach, exponentials are only needed for matrices of small size that depend only on the order of derivatives and number of discretization points, independent of the size of spatial dimensions. This method is particularly advantageous for high dimensional systems, and it can be easily incorporated with IIF to preserve the excellent stability of IIF. Implementation and direct simulations of the array-representation compact IIF (AcIIF) on systems, such as Fokker-Planck equations in three and four dimensions and chemical master equations, in addition to reaction-diffusion equations, show efficiency, accuracy, and robustness of the new method. Such array-presentation based on methods may have broad applications for simulating other complex systems involving high-dimensional data.
منابع مشابه
Karstic water exploration using the Schlumberger VES and dipole–dipole resistivity profiling surveys in the Tepal area, west of Shahrood, Iran
The need for clean groundwater resources to have sustainable development in a country is undoubted. Due to the importance and high quality of karstic waters in supplying water in Iran especially in Shahrood city, it is attempted in this research work to recognize and explore karstic waters in southwest of Tepal area, Shahrood. For this purpose, integration of the results obtained from the metho...
متن کامل2D linear array device as a quality assurance tool in brachytherapy applications
Background: External beam radiotherapy and brachytherapy plays a vital role in the management of cancer cervix. High dose rate brachytherapy is being presently used worldwide for the brachytherapy applications. At present, 2-Dimensional linear array detectors are the most common QA tool used for pretreatment patient specific quality assurance in external beam radiotherapy alon...
متن کاملDirect method for solving nonlinear two-dimensional Volterra-Fredholm integro-differential equations by block-pulse functions
In this paper, an effective numerical method is introduced for the treatment of nonlinear two-dimensional Volterra-Fredholm integro-differential equations. Here, we use the so-called two-dimensional block-pulse functions.First, the two-dimensional block-pulse operational matrix of integration and differentiation has been presented. Then, by using this matrices, the nonlinear two-dimensional Vol...
متن کاملDirect method for solving nonlinear two-dimensional Volterra-Fredholm integro-differential equations by block-pulse functions
In this paper, an effective numerical method is introduced for the treatment of nonlinear two-dimensional Volterra-Fredholm integro-differential equations. Here, we use the so-called two-dimensional block-pulse functions.First, the two-dimensional block-pulse operational matrix of integration and differentiation has been presented. Then, by using this matrices, the nonlinear two-dimensional Vol...
متن کاملOptimized Design of Nanohole Array-Based Plasmonic Color Filters Integrating Genetic Algorithm with FDTD Solutions
Recently, significant interest has been attracted by the potential use of aluminum nanostructures as plasmonic color filters to be great alternatives to the commercial color filters based on dye films or pigments. These color filters offer potential applications in LCDs, LEDs, color printing, CMOS image sensors, and multispectral imaging. However, engineering the optical characteristics of thes...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید
ثبت ناماگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید
ورودعنوان ژورنال:
- Journal of computational physics
دوره 258 شماره
صفحات -
تاریخ انتشار 2014